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On the anomalous t-quark charge asymmetry and 
noncontractibility of the physical space 

Davor Palle 



Heavy flavour production at hadron colliders represents a very promis- 
ing field to test perturbative QCD. The integrated forward-backward asym- 
metry of the top-antitop quark production is particularly sensitive to any 
deviation from the standard QCD calculations. The two Tevatron collab- 
orations, CDF and DO, reported a much larger t-quark charge asymmetry 
than predicted by the theory. We show that the QCD in noncontractible 
space, where the minimal distance is fixed by weak interactions, enhances 
the asymmetry by more than a factor of 3 (5) at the parton level in leading 
order of the coupling for the Tevatron (LHC) center of mass energies. This 
result should not be a surprise since the asymmetry observable directly ex- 
plores the far ultraviolet sector of the spacelike domain of the Minkowski 
spacetime. 



PACS numbers: 12.38.-t; 11.15.Ex 



1. Introduction and motivation 

The missing flavour-mixed light neutrinos, no lepton or baryon number 
violation and the absence of any candidate for a dark matter particle, call for 
a substantial improvement of the Standard Model (SM) of electroweak and 
strong interactions. Namely, the dominance of the baryon over antibaryon 
matter in the Universe suggests that baryon or lepton number be broken in 
particle physics pQ. It is well known, for example, that the broken lepton 
number can induce the breaking of the baryon number [2]. Any alternative 
to the presence of the cold dark matter in the Universe supposes drastic 
changes not only of the General Relativity, but also of the nonrelativistic 
Newtonian theory of gravity [3]. However, the direct or indirect detection 
of the dark matter particle is indispensable. The HESS source J1745-290 
at the centre of our galaxy (for the most recent anaysis see ref. [1]) and 
the hints from the anomalous positron (or antiproton) abundance from the 
PAMELA mission [5] (see, for example, the analysis in ref. [6]) suggest on 
the existence of the very heavy dark matter particle when searching for the 
dark matter annihilation products. 



(2) 
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MOND (with its few relativistic generalizations) represents the alterna- 
tive to the dark matter paradigm (see a discussion in ref. [7]). 

Besides the super symmetric, grand unified or extra dimensional exten- 
sions of the SM, a very conservative alternative to the SM was proposed 
in [8J, called the BY theory, resolving the ultraviolet singularity and the 
SU(2) global anomaly problems. Light and heavy Majorana neutrinos with 
flavour mixing and lepton CP violation could play a crucial role as hot and 
cold dark matter particles in the evolution of the expanding and rotating 
Universe [8j[9]. 

The noncontractible space of the BY theory, as an alternative symmetry- 
breaking mechanism to the Higgs one, introduces into the physical realm a 
new universal Lorentz and gauge invariant constant (UV cutoff in the space- 
like domain of the Minkowski spacetime, see ref. [8]) A = ^ = |-^Miy ~ 
326GeU. 

The enhanced strong coupling at small distances and the absence of the 
asymptotic freedom in QCD are the immediate consequences of noncon- 
tractible space [ID]. We show that electroweak quantum loops with heavy 
t-quark contributing to the CP violating processes of K and B mesons are 
affected by the UV cutoff [IT]. The branching fraction for a rare decay 
B s — > ufi is lower by more than 30% in the BY theory compared with the 
SM owing to the modified short distance part of the amplitude [12] . 

ATLAS and CMS experiments at the LHC reported recently [13] a dis- 
covery of the 125 GeV resonance. It could be the Higgs boson of some 
theory beyond the SM, but it could be also some pseudoscalar or scalar me- 
son with a substantial component of the pseudoscalar or scalar toponium 
[14] . Even the spin 1 boson cannot be excluded as an interpretation of the 
new 125 GeV resonance [15] , 

Anyhow, the Higgs mechanism does not solve the mass problem of par- 
ticles. Eventually, the solutions of the coupled system of nonlinear integral 
Dyson-Schwinger equations of the UV nonsingular BY theory could resolve 
the mass problem of elementary particles [TDJ . 

In this paper, we study the implication of the UV cutoff to the lead- 
ing QCD contribution for the forward-backward asymmetry in the top- 
antitop production. The large discrepancy between the theory and the 
experiment for this asymmetry observable is reported by the Tevatron col- 
laborations CDF and DO [17J. Let us quote the most recent results of the 
CDF collaboration [18] : parton level asymmetry Afb(M^ < 450GeU) : 
Data(±stat ± syst) = 0.084 ± 0.046 ± 0.026 vs. SM expectation = 0.047 ± 
0.014; A FB (M t t > 450GeU) : Data = 0.295 ± 0.058 ± 0.031 vs. SM 
expectation = 0.100 ± 0.030. 

In the next chapter, we present the main ingredients of the calculations 
while providing more details in the Appendix. Results and Conclusions are 
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given in the last chapter. 



2. Charge asymmetry at the parton level 

Almost invariably, various asymmetry observables of the electroweak or 
strong interactions are very sensitive to the details of the underlying pro- 
cesses. It appears that the t-quark pair charge asymmetry can test QCD 
loop corrections [19]. We shall study the dominant quark-antiquark anni- 
hilation channel whose structure equals the electron-positron annihilation 
amplitude modulo coupling and gauge group constant factors [201 ED H2] . 

Let us define the asymmetric part of the differential cross sections [19| 



dojt _ , gg V _ 1 { do(qq -> QX) da(qq ->■ QX) ] 
\ a A ) = 



d cos 6 A 2 d cos 9 d cos 9 

Born cross section (symmetric part of the quark-antiquark annihilation 
to leading order o? s ) is given by [20]. [T9] : 



da(qq -4 QQ; Born) _ 2 T F C F irf3 n 2 2 
d cos 9 JSI C Is 

1 4 / n rn n 

T F =~,C F = -,N c = 3,P=y/l- 4m 2 , m 2 = -2 

2 3 s 

s = £ c 2 m ,c = /3cos0, Z(p(g),p(Q)) = 0. 

The asymmetric part to the leading order consists of the virtual, soft 
and hard gluon emmission differential cross sections [191 \21\ [22] : 



J(/7) o cos 9o\l 1 ok 
(I) 0<k< hi, -1 < cos6> 7 < 1, 
(II) h < k < k 2 , gi(k,E th ) < cos0 7 < g 2 {k,E th ). (1) 

The equations for the virtual, hard and soft gluon radiation in the ap- 
pendix of ref. [19] are obtained from the equations in |21|, [22] in the limit 
of the vanishing mass of incoming fermions. 

The QCD in noncontractible space differs from the standard QCD when 
quantum loops are evaluated with the cutoff in the spacelike domain. Thus, 
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one can find two possible sources of deviation from the standard QCD cal- 
culation for the asymmetry function A°°(cosO) = a' A / 'o' Born : (1) calculation 
of the running coupling oc^ (see ref. [ID]), (2) box diagram contribution to 
the virtual correction [T9 | |2T | |22] : 

(°a)' = (° a)' (virtual, a^) + (a^)' (soft, a^) + (a\)' (dif ference, a^) 



+ (<4)'(hard,t4) = (-|t) 3 (<7a + (a^)' (virtual , o£) - (a^)' (virtual, a j), 



( a Born)' — (~^) 2 ( a B°orn)' ■• 



A «)' ? 
( a Born)' 

(2) 

(a\)' (virtual, a^) — (a^)' (virtual, a^) 

( a Born)' 

A denotes quantity in the BY theory, oo denotes quantity in the SM. 

We mean that (a^)' (virtual, a£) is evaluated with a A coupling, etc. The 
calculation of the strong interaction running coupling in noncontractible 
space was performed in the momentum subtraction renormalization scheme 
to one loop order in ref. [10]. Hard and soft gluon radiations do not contain 
loop diagrams to leading order. 

Our main task should be a reevaluation of the interference term in the 
cross section containing the box diagram in the virtual correction term. To 
accomplish this in noncontractible space, we have to reduce the amplitude 
into pieces that are manifestly translationally and Lorentz invariant. 

We render light quark masses nonvanishing as a regulator of the collinear 
singularity that is canceled away in the asymmetric cross sections. Infrared 
singularity is controlled by the regulator gluon mass and is canceled away 
in both A°° and SA£ asymmetry parameters. The virtual corrections can 
be represented with the following expression [22J 



,4 A fcosi 



A°° + 5 At) + 5 A 



box > 



SAL 



-A c 



°- n box 



G 



(0->7T-0)], (3) 

sJl-Amf/s. 

Definitions are given in the Appendix, as well as the procedure how to 
evaluate the integrals in noncontractible space to maintain translational and 
Lorentz invariance. 

Now we can compare the t-quark charge asymmetries to the leading one 
loop order in the standard QCD and the QCD in noncontractible space. 
The numerics and discussion can be found in the last chapter. 



da A (virtual) _ 3 d 2 abc f3 t ^ 

d* - 4 ° B 

a abc - y> Pt 



3. Results and conclusions 

The difference between the t-quark charge asymmetries of the standard 
QCD and the QCD in noncontractible space lies in the additional two terms 
of Eq.(2) 5A% and 5A£ ox . The first additional term 5A% can be evaluated 
using Table 1 derived from the formulae for in ref. [TU]. This correction 
can enhance the SM asymmetry by up to 47% for the largest parton E cm = 
lATeV. The strong coupling o^(/x) is frozen at \i ~ 0.5TeV. 

This is not enough to explain the asymmetry observed at the Tevatron 
|17j . Fortunately, the second additional term 8A^ ox provides the necessary 
enhancement (see Figure 1 and Table 2). 

We define the integrated charge asymmetry parameter as [19] 



. L a' A dcos0 . A . 
Aint = i • (4) 

Jo a 'Born dcOSe 



One can conclude that the charge asymmetries at the parton level are 
enhanced in the BY theory by more than a factor of 3 (5) for Tevatron 
(LHC) center of mass energies. It is evident from Tables 1 and 2 that the 
deviation from the SM is larger for higher E cm and the virtual correction 
(box diagram) SA^ ox dominates over the strong coupling correction SA^. 
It means that the box diagram explores the deep spacelike domain of the 
Minkowski spacetime to which the asymmetry observable is very sensitive 
and, in addition, there is no new negative compensation of the real hard and 
soft contributions (no quantum loops to this order of perturbation) except 
the new factor. 
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Table 1. Running strong couplings at the scale \i = E cm /2 assuming m u = 
2.5MeV, m d = b.OMeV, m s = WOMeV, m c = 1.6GeV, m b = 4.8GeV, 
ra t = 172GeV and a s (n = M z ) = 0.12. 



E cm (TeV) 


0.4 


1.96 


8 


14 




0.1077 


0.08985 


0.07886 


0.0749 




0.1104 


0.110 


0.110 


0.110 




0.0248 


0.225 


0.397 


0.471 



Table 2. Integrated t-quark charge asymmetries for parton E cm evaluated with 
E th = 0.9 x E cm /2 and m t = \12GeV . 



E cm {TeV) 


0.4 


t 1.96 


8 


14 


4°° 


0.0740 


0.1774 


0.1519 


0.1449 


A A 


0.0939 


0.661 


0.805 


0.874 


AA / 400 


1.27 


3.73 


5.30 


6.03 
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To find charge asymmetry for hadrons, one has to convolve parton cross 
sections with parton distributions. It is necessary to solve DGLAP and 
BFKL equations in noncontractible space. This work remains for the future. 
It is very unlikely that higher orders of perturbation in the strong coupling 
or new parton distributions can remove large deviation of the asymmetry 
from the standard QCD found at the parton level. If the LHC confirms the 
Tevatron results, it will be necessary to investigate the issue to higher per- 
turbative order to reach higher accuracy, because to date, it is the largest 
discrepancy observed between the standard QCD and the experiment. 
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Appendix 

Since the details for the QCD running coupling evaluations can be found 
in ref. [10] and the equations for the SM asymmetries in ref. [19], in the 
Appendix, we outline the equations for the virtual corrections in the SM 
and the BY theory using notations of refs. [21| 122] . 

Let us define the energy unit E = E cm /2 and the dimensionless mass of 
the light quark by rn u = m u /E and the t-quark by nit = Tnt/E |22j . With 
previously defined c = fit cos 6 the coefficients Wj in the sum X)J=i w jlj °f 
Eq.(3) are as follows [22] : 



ioi = 1 + c 2 - 2c 3 + (1 - 2c)(m 2 + m 2 ), w 2 = 2c(l - c) - m 2 - cm 2 , 



ws = — 1 — c, wq = 1, W7 = 1 — c. 

We need further definitions to describe the process q(p+) + q(j>-) — > 
t(q + ) +i(q-) and its amplitude [21] 122]: 



W3 = 2c(l — c) — m t — cm u , W4 = 2 — c + c + m t + m: 




(A) = k 2 - 2k ■ A - P 2 + is, (Q) = k 2 - 2k ■ Q - P 2 + is, 
{±) = k 2 ±2k-P + P 2 -m 2 gluon + ie. 

Now we can define dimensionless integrals Ij of eq.(3): 
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h = E 4 J (1), h = E 2 J (k ■ A), I 3 = E 2 J (k ■ Q), h = E 2 J (k 2 ), 

h = J ((k ■ P) 2 ), h = J {{k • A) 2 + (k ■ Q) 2 ), I 7 = J ((k ■ A)(k ■ Q)). 
These integrals can be evaluated by the integrals from ref. |21j 



\t. t . T 1 - f rfih [lijjfi^d 
L ' M ' VJ ~ J d " (A)(Q)(+)(- 



(A)(Q)(+)(- 

that are expressed in terms of nine functions: F, G, Fa, Fq, C?a, Gq,Hp, 

Ha, HQ- 
Let us represent seven integrals Ij of ref. [22J in terms of functions from 
EU: 



h = -^(A 2 J A + A-QJ Q )E 2 1 
h = ^(Q 2 Jq + A-QJ a )E 2 , 

/ 4 = - 4 

7T 

h = ^(KqP 2 + Kp(P'Y), 

4 
TT 

4 

TT 

where J a, Jq, Ko, Kp, K A , Kq and Kx functions are defined in 
terms of nine functions F, . . . , Hq |21| . 

We use standard definitions for the scalar two, three and four point 
functions |23|: 



2 (4K + P 2 K P + A 2 £T A + Q 2 K Q + 2 A • QK X )E 2 , (5) 

4 ^2 r ^ /d2n2\ 



J e = 4(^o(A 2 + g 2 ) + E:A((A 2 ) 2 + (A-Q) 2 ) + ^Q((A.Q) 2 + (g 2 ) 2 ) 

TT* 

+ 2F/ X A-Q(A 2 + Q 2 )), 
h = ^(A • QF/ + A • QA 2 K A + A • QQ 2 K Q + K X (Q 2 A 2 + (A • Q) 2 )), 



B (p;mi,m 2 ) = (wr 2 ) -1 / d 4 k[k 2 - mf + ie] 1 [(k + p) 2 - m\ + le] \ 

x[(^ + P2) 2 -m-| + «e] _1 , 
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D (p 1 ,p2,P3;rno,m 1 ,m 2 ,m 3 ) = (in 2 ) 1 J d 4 k[k 2 - ml + le] 1 
x [(*: + Pi? -m\ + if]' 1 [(k + p 2 ) 2 -m 2 2 + ie]~ 1 [(A; + p 3 ) 2 - m\ + ie] _1 . 

In ref. [21] expressions for all nine functions F,...,Hq in the standard 
QCD can be found. The same functions have to be expressed by the previous 
scalar two, three and four point Green functions in noncontractible space 
in order to properly restore translational invariance |10l [TT] IT2] . 

Functions G,F a ,Fq have already a suitable form of the three point 
functions |21|: 



G = J d 4 fe(A)- 1 (Q)- 1 (+)- 1 , F A = j d 4 k(A)-\+)-\-)-\ 

f q = J <£k{Q)-\+y\-y\ 

Note that all expressions in ref. [21] are derived under the assumption 

of 

fngiuon = A <C m u ,m t , E cm . From their definitions, Ga and Gq can be 
expressed as: 



% G Q = 7fl^ F Q + ^S[^o(-2P; A, A) - UB (-Q - P; A, m t )}, 
Pt s Pt 
1 2vr 2 

%G A = -^%F A + _[SRB (-2P; A, A) - UB (-A - P; A, m,)]. (6) 

Pu s Pu 

For functions F, Hp, Ha, Hq, we derive the equations that allow to put 
these functions in the alternative form expressed only through scalar n-point 
integrals. 

The linear system for the F function looks as 



pirn + pi ■ P2m + pi ■ P3m = Ri, 
pi ■ P2V1 + plm + P2 ■ P3m = R2, 
pi • P3V1 + P2 ■ P3V2 + plm = R3, (7) 



R i = ^[KC (p2,P3;"io,m2,m 3 ) - KC (p 2 ~Pi,P3 - Pi; mi, m 2 ,m 3 ) 
-(Pi - mf + rnl)^tD {p 1 ,p 2 ,P3-,rno,mi,rn2,rn3)], 
R 2 = -^[^C Q (pi,p3;mo,m 1 ,m 3 ) - $tC (p 2 ~Pi,P3 ~ Pi\ mi, m 2 , m 3 ) 
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-ipl - ml + mo)5ft-Do(Pi,P2 ) P3;mo,mi,m2,m 3 )], 

^3 = -[$tC (pi,p2;mo,m 1 ,m 2 ) - $tC (p 2 ~Pi,P3 - Pi;m 1 ,m 2 ,m 3 ) 

-(pi ~ m i + ml)^D (pi,p2,P3;mo,m 1 ,m 2 ,m 3 )], 
Pl = 2P, p 2 = P - A, p 3 = P - Q, m = mi = X, m 2 = m u , m 3 = m t 

^%F = -QF Q + 2vr 2 (A 2 r ?2 + A • Qrj 3 ). 

Similarly, we derive the linear system for H functions 



P1P1 +Pi -V2P2 = Mi, 

Pi -P2Pi + P 2 P2 = M 2 , (8) 

Mt = ^B {p 2 ;X,m 2 )-^B {p2-pi;m 1 ,m 2 ) + (-X 2 + m\-p\) 

x ^C (p 1 ,p 2 ;X,m 1 ,m 2 )], 

M 2 = -[m (p 1 ;\,m 1 )-m (p 2 -p 1 ;m 1 ,m 2 ) + (-X 2 + m 2 2 -p 2 2 ) 

x i RC (p 1 ,p 2 ;X,m 1 ,m 2 )], 
Pi = P ~ A, p 2 = P - Q, m 1 = m u , m 2 = m t 

=> ZHp = QG + ir 2 ( Pl + p 2 ), ZH A = -7r 2 pi, 9F Q = -7r 2 p 2 . 

The validity of new forms for F, Hp, and Hq is also checked numer- 
ically. 

The virtual corrections can be evaluated by eq. (A.l) of ref. [22J or by 
eq. (12) of ref. |2I]. 

We are now prepared for the crucial step to calculate virtual corrections 
in noncontractible space defining scalar n-point integrals in noncontractible 
space. Bq function is outlined in refs. [10j [HI 02] . The similar procedure 
should be applied to the three point function: 



fiC °° = 5RC A + 6Cfr(symm), (9) 
^(symm) = -[SCq(pi,p 2 ; m , mi, m 2 ) + 5Cq(~pi,p 2 - py, mi, m , m 2 ) 
+5Cq(-p 2 ,P! - p 2 ;m 2 ,m ,mi)}, 

rl/A r+1 r+1 

SC (px,p 2 ;mo,mi,m 2 ) = tt~ 2 / dww~ 5 / dx\J 1 — x 2 / 

J «/-! J-l 
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i-2-k 

x / d(p[-k 2 -mgj-^-fc 2 + 2{k-pi)+ P l-m\}- 1 
Jo 

x [-A: 2 + 2{k ■ p 2 ) +pl- rr%]-\k = w' 1 ), 
where (k ■ p{) = tkx(pi)° — k ■ pi, 

k = ky/ 1 — x 2 {\J 1 — y 2 cos <f>, \Jl — y 2 sin (/), y). 

All the imaginary parts of the subintegral function in 5Cq are erased by 
integration as odd functions in variable x. 

The same decomposition is possible for the four point function, although 
with four terms necessary for symmetrization in SDq. 

Multidimensional numerical integrations in virtual and real gluon radia- 
tions are performed by Suave routine from CUBA library [23] to the relative 
accuracy of 0(1O -4 ) with up to 50 million of sampling points per integral. 
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